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ABSTRACT.—It has been suggested that archaeologically recovered cottontail mandibles 
may be identified to species by examining the relationship between mandibular depth and 
alveolar length even though most of the cottontail skeleton is not diagnostic to species. 
Study of cottontail mandibles recovered bv the Dolores Archaeological Program from Ana- 

w 

uzi sites in southwestern Colorado has supported this suggestion but indicated that reliable 
species identification requires a more thorough analysis of traits. Further investigations 
should be directed to the topic of subspecific variation as well. 


INTRODUCTION 


Although the skeletal remains of cottontail rabbits commonly are recovered from 

archaeological sites (e.g. Bertram and Draper 1983, Binford et al. 1983, Cordell 19/7), 
Geological separation of these cottontail species seldom is possible. Findley et al. 
19 * 5 . 83 - 86 ) have noted that the taxonomic relationships between the eastern cotton 
^ Sylvilagus tlondanus), the mountain or NuttalTs cottontail (Sylvildgus nuttallu), and 
^ desert cottontail (Sylvilagus audubonii) require further clarification, at least in New 

\|ij X iC 0 v ’ i _] rVior. 


Nevertheless, these authors have sugge 


mandibular 


-uicac aumors nave suggcaicu ui#i -- 

can used to separate the desert cottontail from the other two species e\en in 
** ac °’°gical and fossil materials. These characters, particularly mandibular measure- 

TTiPrL _ 1 . m __*«*•««* in 



^ m e the species represented (Akins 1984; Anderson 1980; Harris 196 . > 

m 1979 )- This article indicates that the identification of cotton tails is not neces- 
as straightforward a matter as suggested by Findley and his coworkers. 

DoW C k° nes °f cottontail rabbits occur more frequently in sites excavated > 

For tv.* ‘^ rc ^ ae °logical Program 1 than those of any other single taxon (Neusius 
ttiirc S rcason a ^°ne it was considered desirable to determine the cottontai sp 

nte d. However, because there was reason to expect both desert an ut 
to be present and because the habitat preferences of these species wer 

'Armstrong 1972:82, 85; Bissell and Dillon 1982:7), further study o e 


Dol 
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Anasazi 


Neusius in press). Thus, we used the characters 


ti-,. ;V 5 UJ assi S n cottontail mandibles recovered at Dolores to specie . . 

^ 1Ca hon ii * vaiaate plots less than satisfactory because the accuracy o spec ^ or 

our purr n0t ^nt. While the use of discriminant analysis was more satisfac Y 

(Sblej° SeS ’ a PP r °ach suggested that the proper identification o cot on 

Itl0re CorT1 plex task than anticipated and deserving of further stu v- 

SPECIES DETERMINATIONS USING BIVARIATE PLOIS 

Cr? ,0 Findle V et al. (1975:84-85) the desert cottontaU has a deeper j aw an da 

throw than either tbe M.,troii’c the extern cottontail. The charac 
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to measure this are the alveolar length and the depth of the lower jaw (Fig. 1 ). Alveolar 

length is defined as the distance between the anterior margin of the third premolar and 

the posterior margin of the third molar. However, these authors suggest that the distance 

between the anterior border of P 3 and the posterior margin of Ml will suffice in studies 

of archaeological and fossil material. Mandibular depth is defined as the distance between 

the anterior border ot the fourth premolar and the ventral border of the mandible at right 
angles to the toothrow. 


Dental characteristics also may be useful in identifying cottontail species according 
to Findley and his coworkers. Ihe borders of the desert cottontail’s posterior externa, 
reentrant angle on premolar 3 are strongly crenulated, while those of the eastern cotton¬ 
tail are less crenulated and those of Nuttall’s cottontail are smooth (Findlev et al. 1973: 

85-86).3 

In the Dolores study osteological measurements were made on 302 cottontail man¬ 
dibles from 27 archaeological sites located in the Dolores project area. Nineteen of these 

n r\ IT rk + 1 « • C _ • *1 » « 1 ■ « - _ 


analysis. 


an 


Ihe measurements taken followed those recommended by Findley and his co- 

workeis. As shown in Figure 1, distance A is the depth of the jaw, measured with a ver¬ 
nier calipers. In order to standardize this measurement, the mandible was placed on 
graph paper with the cheektooth row aligned vertically and the labial side up. The 
stationary part of the calipers then was placed on the most anterior portion of the alveo¬ 
lar notch between P 4 and P 3 and the sliding portion was brought to the basal border of 


-- - iu li ic aligned cneeKtootn row. However, uie ucvuupn.-.. 

of the ridge ot bone at the anterior alveolar notch is variable and the area may be eroded 
in archaeological specimens. As with all our measurements we attempted to be as consis¬ 
tent as possible in locating this point. Distance B is the length of the toothrow. The 
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the mandible. „„ lunc 

alveolar length (Pj-Mj). 


the depth of 
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(P3-M3) and distance c n 


the cottontail mandible. Distance A is 
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measurement was taken with the calipers from the anterior alveolar notch of the third 
:'emolar to the posterior alveolar notch of the third molar. Distance C is the modified 
alveolar length recommended for archaeological materials. The calipers were set at the 
r.terior alveolar notch of the third premolar and expanded to the posterior notch of 
the first molar in order to obtain this measurement. In our study, sample size was in¬ 
creased by 27 when mandibles for which only distances A and C could be measured 
were added. 

The crenulation of the enamel border of the third premolar also was examined. 

Using a 10 power magnifying lens, a subjective evaluation was made on a scale of 0 to 

3 with 3 assigned to very-, 2 to somewhat-, 1 to slightly crenulated, and 0 to uncrenulated 
specimens. 

In order to classify these specimens each mandible for which distance B was present 
was plotted on a scattergram (Fig. 2). Mandibular depth was measured along the x-axis 
-id alveolar length (P 3 -M 4 ) was measured on the y-axis. An approximation of the line 
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separating desert and NuttalTs or eastern cottontail provided by Findley et al. (19" 
Fig. 35) was then placed on the scattergram. Mandibles falling above the line were identi 

fied as Nut tall s cottontail while those falling below the line were identified as desert 
cottontail. 


A second scattergram was constructed for the relationship between mandibular 

depth and the shorter measure of alveolar length (P 3 -M 1 ). In this scattergram the species 

assignments made on the basis of the first scattergram were used to assess the probate 

species of the mandibles for which distance B was lacking. In those cases in which .1 

remained unclear, the crenulation of the enamel border was used as the arbiter of pr 

able species. However, we considered these identifications to be unreliable compared tt 
those made obvious by the bivariate plot. 

Ut then used the species determinations to test the expectation of localized pre 

curement of cottontails among the Dolores Anasazi by examining distribution of species 

across space and through time at Dolores archaeological sites. This has been reported! 
Dint and Neusius (in press). 


The 


Ho\\e\er, we remained dissatisfied with the method of species determination 

tt • • ^ - A. - 1 tt • A A _ _ 



WVIU jyv.uv.3 VVA, I L piLJCUL 111 Lilt LJKJ. 

uti lze< ^ the Anasazi. They did not indicate how reliable the species i 


i 


1 iu 



were. 


In instances in which our subjective evaluation of crenulation of the enamel bord 
was the primary indicator, we considered the identifications unreliable, but more subtle 

_ ® . • £ w 

'aria ion in accuracy was not obvious. Because correct identification was important to 

analysis of the temporal and spatial distributions of these species, we wanted to assess 
our accuracy more thoroughly!, 


SPECIES DETERMINATIONS USING DISCRIMINANT ANALYSIS 


us, we made a second attempt to determine species using mandibular measure* 
nts. We decided to use discriminant analysis on a group of modem cottontails \\hose 

species were known _ i___ ,. . . „ . _nr^his* 


toric unknowns. 


Uim ultIi employ me aiscriminant runction to ciassuy r*' 
Discriminant analysis is an appropriate technique for such classification 

-u), and provides probabilities for the group assignment of each case. 

, C anc ^ Wildlife Service in Fort Collins, Colorado provided a sample 0 . 

nown cottontail mandibles from their skeletal collections. Distances A and B tor 4- 


Sendee 


b already had been measured, and we wen 
same indhiduals. Twenty-five of the individuals 

_ M ■ A A _ 


distance C 


while 


t , , -- uuuuuunn ualleyi ai uic unit '-“i— 

on f . . Cen as Sylvilagus nuttallii grangeri. We also had measureirun 

T ! V Vn UOnal ind ividuals from die collections at Mesa Verde National Park and For. 
t^cwis Co eap urki^u u.j . these 


Ipu-ic n 11 . luucuiuiis ai mesa verue 

had h ° • j^ C V ^. ^ a d F> ee n identified from soft tissue characteristics. One o 
Hall iqs! ^ntified as Sylvilagus audubonii warreni (Sylvilagus audubonn ccdrophilus? 

simnlv • a ? ^ ne was considered Sylvilagus nuttallii pinetis and three wf 

whirh aSSlgne , d to S y lv Vagus audubonii. This gave us forty-six known individuals«'» * 

WOrk - Alt hough this sample was not overly large, it was larger than that ip- 

th< 


parentlv used by Findley et al. (N=35). 


Two discriminate function 


an 


these individuals uSin!f , 

— - al. 1975:434 


to the m ^ ese P er f°rmed a stepwise selection of three variables correspo ^ 

lated witrdilt'ancrB 10 *-' 1 *'’' 31 !" 5 A ' 8 and C Because distance C WaS ^letyon' 

fcncUon » P " Cent ° f tha ">odem individuals were correctly classified us,ng * 

uucu °n. because fhoca.- 1_;_. . .... - ... JwMncC c 


mandibular fragments 


using 


tenths percent o~f t h V “ tc<15uremen « tor distances A and C. Ninety 

e nowns were correctly classified using this function. 


-three 


and 
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Scores on the second discriminant function were then used to assign species to the 
L3 adult Dolores mandibles originally classified using bivariate plots. Figure 3 is a 
-..stogram displaying the distribution of cases with respect to this second discriminant 
function. This procedure provided us with probabilities for group membership. 
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CANONICAL DISCRIMINANT FUNCTION 1 
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set jailing to the left of 0 were 

**assigned to desert cottontail. 
' f discriminant scores 
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Probabilities for these assignments vary according ti 
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DISCUSSION 


ait u&i|jui tail l iu Aw 

u *«i desert cottontails tend to have a deeper and longer mandible tnan .\uium a 
bas d mai!S ^ suggested by Findley et al. (1975), but also indicates that classification 
>lk°L y ° n this rela tionship is less than ideal. If we consider a probability of group 
could * J* .°^ ad equate for identification, 35.8% of the prehistoric man i 
beidentifi d ^ en ^^ e< ^‘ ^ a more stringent requirement of .95 is emploced 45.7/o cann 

’n (J le [' Inore ’ t b e probabilities for those mandibles assigned to Nuttall s 

ie °' ver than those for mandibles assigned to desert cottontail. 0 \erhalf Io 

^ h ibleS assi § ned to Nuttall’s cottontail have probabilities of less than .901 while 
Waller n C i° r ° n ^ of those assigned to desert cottontail. This may resu t ro 

this sperj(» mber ° f known Nuttall’s cottontails or it may indicate greater variability wit m 


First, our analysis con- 


rl . 

Vi 
*2 * 


i rehi'? tereSt * n ' 8 ’ l ^ lat desef t cottontails apparently represent such a high proportion 
«itif auna , 0nc man dibles. We expected Nuttall’s cottontail to predominate in e 
^ (Bve . assem blages both because much of the project area was foreste prc “ 

National' 1 -° ? ress ’ ^ cte rsen, in press) and because much of it was above e 

fre obsen- l mU for deser t cottontail (Armstrong 1972:85). However, since our su J j 

no att l ° nS the project area today did not fit with this expectati . 

b’rttil better if^ >t . We *§ht the discriminant analysis in favor of one or the ot e P 
tiip n_i ls crimination .i,ic r.r»ccihle. we can only specu a 


tics. 


important to understanding 
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T he discriminant analysis also results in some mandibles being assigned to a differ-n 
species than was suggested by the use of bivariate plots. Thirty-nine (10.2%) of the pr 
historic mandibles were assigned to a different species by the discriminant analysis. I 
probabilities of group membership do tend to be low for this group. Within it 74.4’oot 
the cases are assigned to a species with probabilities of less than .90. Nevertheless, it is 
clear that the modern population we obtained differs from that used by Findley and his 
coworkers with respect to these two characters. Sample size could be a factor here, bi' 
it is also possible that subspecific variation is involved. Most of the modem desert cotto: 
tails measured by us were assigned to S.a. bailey i, a subspecies which does not occur :n 
New Mexico (Hall 1981 :Map 225). Similarly most of the modem Nuttall’s cottontails 
were assigned to S.n. grangeri which also does not occur in New Mexico (Hall 198F 

Map 223). 

Finally, the use of discriminant analysis in zooarchaeological analysis may be pre¬ 
ferable to that of bivariate plots even if such a small number of variables is involved. 
I his is because it provides probability assessments for the species assignment ot each 

mandible. However, the discriminant analysis also indicates that the relationships be- 

# 

tween these species is not adequately defined by these characters alone. Separation 
could be improved by obtaining a larger sample of known specimens, adding more char* 
acters to the analysis^ and investigating the possibility of subspecific variation. 

We intend this article as a caution to zooarchaeologists attempting to identify cotton¬ 
tail species. Certainly it is naive to apply the plot provided by Findley et al. toavarkt* 
ol archaeological situations without further research on the variability in these characters 
between species. However, we do not mean to suggest that the separation of cottont $ 
not be attempted by zooarchaeologists. A number of important behavioral and environ 
mental questions may be addressed through detailed study of cottontail remains, 
think that this is an instance in which further investigation of taxonomic and ciistrili 
tional relationships by zooarchaeologists is required in order to do justice to zooarcha 
ological data bases. We hope that this summary of our analyses serves as a stimulus 
further study of this topic. 
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NOTES 


Since 1978 the Dolores Archaeological Program has been under contract with the l\> 
Department of Interior, Bureau of Reclamation (Contract No. 8-07-40-S0562) to mil 
gate the effects of the construction of the McPhee Dam and Reservoir in Montezuma 
County, Colorado. Faunal remains have been recovered from approximately 150 archj 
ological sites the majority of which belong to the Anasazi cultural tradition and date 

between A.D. 600 and 950. 


" The modem range of eastern cottontail is far removed from the Dolores area (Hall 1981- 
Map 223). Ihus, we did not anticipate the presence of this species in the area during th? 
period of Anasazi occupation. 

Crcnulation oi the enamel of all cheek teeth has been examined by Orr (1940) for 
California cottontails and may be useful for Colorado cottontails as w r ell (R. Findley pers. 

communication). 


4 


I he authors are aware that more variables are normally used in discriminant analysis. 
However, Klecka (1980:15) specifically indicates that we have not violated any basi 
assumptions of this technique. We used this technique because it allowed us to assign 
group membership or species with a probability, the exact problem w r ith the bivariate 
plots which we sought to remedy. 


Hulbert (19/2) uses dental characters which might be of use. 



